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Statement  of  the  Problem 


A  three-year  reeearch  program  to  investigate  the  performance  and 
potential  of  the  laser  time-of-flight  (LTV)  valocimeter  with  photon-burst 
correlation  detection  for  remote  sensing  of  atmospheric  crosswinds. 

Summary  of  Results 

i.  The  scaling  law  for  the  performance  evaluation  of  crosswind 
measurements  with  visible  velocimeters  has  been  established 
theoretically  and  confirmed  with  laboratory  and  field  experi¬ 
ments. 

11.  The  effects  of  atmospheric  turbulence  on  laser  beam  separation 
and  beam  radius  have  been  assessed  experimentally  with  a  diode 
array  and  the  results  interpreted. 

ill.  The  signal  particle  arrival  rate  for  crosswind  measurements  at 
a  range  of  500  m  has  been  predicted  from  the  scaling  law  and 
the  results  agreed  with  the  measured  rates  within  the  experi¬ 
mental  errors. 

iv.  With  the  successful  measurements  at  a  500  m  range,  reported  for 
the  first  time,  it  la  suggested  that  with  proper  method  for 
signal  analysis,  routine  operation  of  visible  laser  valocimatry 
at  long  ranges,  up  to  1000  m,  should  be  feasible. 
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I.  Introduction 

It  la  well  known  that  the  frequency  of  scattered  light  from  a  moving 
particle  la  Doppler  shifted  by  an  amount  proportional  to  its  velocity. 

This  concept  was  first  used  in  laser  light  scattering  to  meaHiire  flow 
velocities  in  1964  by  Yah  and  Cummins  [1].  Since  that  time,  variations 

t 

of  their  technique  have  been  suggested  theoretically  and  tested  experi¬ 
mentally  in  the  laboratory,  as  well  as  in  the  atmosphere,  with  consid¬ 
erable  success.  A  very  recent  review  by  Danielsson  and  Pike  [2]  has 
provided  an  excellent  general  account  on  the  history  of  laser  valodmetry, 
the  comparison  between  coherent  and  incoherent  systems,  and  the  feasi¬ 
bility  for  long-range  atmospheric  wind  sensing.  We  concur  with  the  main 
conclusions  of  this  review  that  because  of  the  large  coherent  detection 
area  at  10.6  urn  and  its  high  laser  efficiency  and  power,  the  coherent 
C02  Ildar  system  is  most  aultable  for  long-range  radial  wind  velocity 
measurements,  but  that  when  higher  spatial  resolution  is  required,  laser 
valodmetry  using  visible  wavelengths  with  photon  correlation  analysis 
still  offers  the  best  results  fot  atmospheric  crosswind  measurements. 

In  addition  to  aircraft  applications,  intensive  studies  for  a  Global 
wind  measuring  satellite  system  (WINDSAT)  based  on  CO^,  coherent  Ildar 
technology  are  underway  [3]  and  the  routine  operation  of  its  ground- 
based  prototype  system  has  been  reported  [4],  On  the  other  hand, 
although  impressive  progress  in  atmospheric  crosswind  measurements  with 
visible  wavelengths  and  photon  correlation  has  been  made  [5],  further 
developments  and  investigations  for  long-range  applications  are  still 


needed . 
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In  this  report,  we  first  review  the  progress  of  atmospheric 
crosswind  measurements  In  the  visible  for  ini  'mediate  ranges  (roughly 
25  m  -  150  m)  with  Che  emphasis  on  photon-burst  correlation  techniques 
which  we  developed  for  Individual  particle  scattering  [6-8].  The  effect 
of  atmospheric  turbulence  on  beam  radius  and  on  performance  evaluation 
and  prediction  for  long-range  applications  from  the  scaling  law  [9,10] 
is  discussed  both  theoretically  and  experimentally.  Finally,  the 
prospect  for  long-range  (roughly  150  m  -  1000  m)  crosswind  measurements 
with  cw  lasers  In  the  visible  wavelengths  using  photon-burst  techniques, 
perhaps  modified  to  handle  the  unavoidable  multiple  particle  scattering 
in  these  ranges  Is  suggested. 

Although  quite  related  to  each  other,  the  optical  system  for  cross- 
wind  measurements  in  the  visible  is  commonly  classified  as  the  differ¬ 
ential  Doppler  (LDV)  and  laaar  transit  or  time-of-f light  (LTV)  systems. 

A  brief  listing  for  successful  atmospheric  crosswind  msasuremsnts 
without  artificial  seadlng  according  to  ths  range  and  year  of  the  work 
rune  as  Bourke  and  Brown  (30  m,  1971  [11]),  Farmer  and  Homkohl  (30  m, 
1973  [12]),  Burtlett  and  She  (60  m,  1976  [6]),  Danielsson  (60  m,  1980 
[13]),  and  Durst  at  al  (105  m,  1980  [14])  with  LDV,  and  Bsrtlstt  and  She 
(100  m,  1977  [7])  and  Lading  at  al  (70  m,  1978  [15])  with  LTV.  We  note 
that  photon  correlation  is  the  key  to  success  of  all  this  work,  The 
concept  of  single-particle  scattering  and  single  photon-burst  correla¬ 
tion,  which  Is  knowingly  or  unknowingly  involved  in  the  crosswind 
measurements  for  intermediate  ranges,  has  been  the  main  theme  of  our 
work  [6-10]  and  will  be  reviewed  here  briefly.  We  also  noticed  thAt 
more  recently,  Durst,  Howe,  and  Richter  [16]  have  reported  a  complete 
LDV  system  with  automated  data  acquisition  suitable  for  routine 
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measurement  of  crosswind  velocity  component  at  Intermediate  range 
(106  m) .  This  system,  whose  data  analysis  is  capable  of  handling 
photon  bursts  resulting  from  multiple-particle  scattering,  as  well  as 
single  particle,  may  be  suitable  for  applications  in  longer  ranges. 
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II.  Single  and  Multiple  Particle  Scattering 

Natural  aerosol  comes  in  different  sizes  and  at  random  time 
Intervals.  Depending  on  spatial  resolution,  incident  laser  power,  and 
detection  system  which  may  be  set  to  partially  eliminate  weaker  back¬ 
ground,  the  signal  for  a  velocity  measurement  consists  of  a  single  or 
multiple  photon  burst  resulting  from  single  or  multiple  particle  scat¬ 
tering  [8].  The  performance  of  L If  and  LTV  depends  on  aerosol  conditions 
and  the  manner  in  which  these  photon  bursts  are  analyzed. 

Perhaps  the  simplest  way  to  appreciate  the  differences  between 
■ingle  and  multiple  particle  scattering  la  to  consider  several  hypo¬ 
thetical  situations  of  light  scattering  from  dust  particles,  as  depicted 
in  Fig.  1.  Four  different  combinations  of  dust  particles  are  shown  to 
travarsa  a  set  of  interference  fringes  on  the  left  and  two  illuminated 
laser  beams  on  tha  right.  Case  (a)  conaists  of  many  amall  particles 
which  scattsr  an  almost  constant  intensity  of  light  n(t),  aay,  in  an 
arbitrary  unit  of  10,  at  all  times.  The  small  intensity  variation 
detectad  in  this  case  is  the  result  of  fluctuations  of  Polason  statistics. 
As  can  ba  seen  in  case  (a) ,  the  simultaneous  scattering  from  many  small 
particles  at  random  positions  yields  no  uaeful  speed  information.  If  a 
single  particle  big  enough  to  yield  a  burst  of  light  sufficient  for  a 
speed  measurement  (we  shall  call  it  a  signal  particle)  traverses  the 
interference  fringes  and/or  the  two  illuminated  baams,  the  detected 
scattering  Intensity  maps  out  the  intensity  profile  of  the  viewing 
volume  and  provides  the  information  necessary  for  the  wind-speed  meas¬ 
urement.  As  shown  in  case  (b)  in  Fig.  1,  the  scattered  intensity  is 
much  less  than  the  detected  background  light  of  case  (a) .  but  itB 
temporal  variation  contains  the  desired  speed  information,  A  physically 
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very  unlikely  case  is  case  (c) ,  in  which  five  signal  particles  are 
aligned  in  parallel  to  the  interference  fringes  and  are  traversing  the 
viewing  volume  together.  This  is  very  desirable  because  the  contribu¬ 
tions  of  these  particles  even  in  the  incoherent  case  add  to  produce  a 
five-fold  Increase  in  the  intensity  variation  for  speed  measurements. 

The  small  fluctuations  in  intensities  in  cases, (b)  and  (c)  are  again 
statistical.  In  case  (d),  20  randomly  distributed  signal  particles  are 
assumed  to  traverse  the  scattering  volume.  The  background  light  is  more 
Intense  then  that  due  to  one  large  particle,  but  the  variation  In 
intensity  which  contains  the  velocity  information  is  proportionally  much 
less  due  to  Interferences  resulting  from  the  random  positions  of  these 
scatterers.  Aa  thasa  particles  traverse  the  two  illuminated  beams,  the 
scattered  Intensity  bursts  are  much  wider,  resulting  in  a  partial  wash¬ 
out  of  speed  information.  It  is  apparent  that  the  signs 1-to-noise  ratio 
resulting  from  the  20  scatterers  (case  d)  is  no  better  than  that  dua  to 
a  single  signal  particle  (case  b).  In  fact,  as  the  number  of  scatterers 
continue  to  increase,  the  visibility  of  the  Interference  fringes  and 
illuminated  beam  profiles,  as  measured  by  scattered  light,  would  be 
washed  out  end  the  desired  speed  information  no  longer  contained  in  the 
detected  light  intensity.  Therefore  the  signal  for  a  speed  measurement, 
when  feasible,  is  best  provided  by  the  scattering  of  a  large  single 
particle. 

Although  multiple-particle  scattering  degrades  slgnal-to-nolse  for 
speed  measurements  as  depicted  in  Fig.  1,  a  process  called  clipping  used 
in  performing  digital  correlation  [8,17],  which  eliminates  the  background 
of  case  (d)  in  Fig.  1,  will  enhance  the  slgnal-to-noise  considerably 
unless  the  number  of  scatterers  is  truly  overabundant,  as  depicted  in 
case  (a)  of  Fig.  1.  Therefore,  a  photon  burst,  resulting  from 


simultaneous  scattering  from  a  few  particles,  when  properly  processed  is 
also  useful  for  speed  measurement.  Such  a  photon  burst  can  be  important 

and  ahould  be  Included  when  long-range  crosswind  measurements  are 
considered. 
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III.  The  Photon  Burst  Correlation  Method 

To  gain  eon*  feeling  about  the  typical  photon  burst  sires 

encountered  in  a  crossvlnd  measurement,  we  consider  the  signal  frem  a 

dust  particle  2  pm  In  diameter  backacattered  from  a  0.1  W,  0.5  cm 

diameter  laser  beam  at  a  distance  of  20  m  from  the  detector.  The  back- 

2 

scattering  differential  cross-section  [18}  Is  approximately  0.2  r 

where  r  is  the  radius  of  the  particle.  Since  the  laser  intensity  is 

0.5  W/cm  ,  the  differential  back- scattered  power  would  be  (0.5)(0.2)(1  x 
-4  2-9 

10  )  ■  1.0  «  10  W/sr.  With  a  20  cm  diameter  receiving  telescope, 

the  detection  solid  angle  is  7.85  x  io“*  sr.  For  a  laser  at  5145  X  and 

a  photomultiplier  tube  with  20  per  cent  quantum  efficiency,  the  detected 

4 

photon  rate  is  calculated  to  be  4  x  10  counts  per  second.  Xf  the 
particle's  transverse  velocity  is  100  cm/s,  it  stays  in  the  0.5  cm  beam 
for  5  ms,  giving  rise  to  a  total  detection  of  200  photon  counts.  This 
is  quite  adequate  for  obtaining  a  single-burst  digital  correlation  to 
determine  wind  speed,  In  practice,  more  than  0.1  W  of  power  could  be 
used  to  enhance  the  desired  single-particle  signal  for  measurements  at 
Intermediate  ranges.  At  long  ranges  (>150  m),  photon  bursts  resulting 
from  multiple-particle  scattering  should  also  be  considered. 

The  photon  burst  signal  is  analyzed  by  a  digital  correlator,  in 
which,  time  is  divided  into  discrete  sample  times,  and  the  number  of 
photon  counts,', t* (i)  in  each  sample  time  is  measured.  The  digital  corre¬ 
lation  function  of  this  signal  for  a  sequence  of  N  samples  is  defined  as 

N 

R(t)  -  R(j t  )  -  I  n(i)n(i-j)  (1) 

i-j 


J  %  V  r.  < i  v 
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where  1  <  J  <  J  and  J  is  the  maximum  number  of  delay  channels  in  the 
correlator.  To  speed  up  the  computation  of  R(t),  a  process  called 
■ingle  clipping  [17]  is  often  used.  In  this  method,  n(i»j)  is  set  to 
either  1  or  0  in  a  shift  register  depending  on  whether  or  not  the 
actual  n(i-J)  Is  greater  than  k,  a  preset  Integer  called  the  clip 
level.  The  time-consuming  multiplication  in  (1)  can  then  be  done  by 
first  ending  each  pulse  of  the  present  sample  n(i)  with  the  shift 
register  contents  and  then  accumulating  the  results  in  the  memory 
stores,  to  form  R(t)  in  reel  time.  When  a  relatively  large  photon 
burst  is  analysed,  n(l)  exceeds  the  clip  level  in  most  samples  during 
the  burst  so  that  n(l-J)  ■  Is  the  single-clipped  correlation  function, 
extrapolated  to  R(0)  at  t  •  0,  gives,  at  least  approximately,  the  total 
number  of  counts  in  the  photon  burst.  The  size  of  the  photon  buret  (In 
counts)  can  than  be  reed  directly  from  the  correlation  function  et 
t  *  0  is  R(0)  minus  the  noise  background. 

The  experimental  arrangement  utilized  in  the  LTV  field  measurements 
is  shown  in  Fig.  2.  The  output  of  a  1.5  W  cw  argon-ion  laser  at  5145  X 
is  sent  to  a  Malvern  model  RF307  beam  splitter.  Using  a  single  telescope, 
the  two  beams  from  tha  beam  splitter  ere  focussed  at  the  sensing  volume. 
The  scattered  light  is  collected  by  a  Celeetron  8-inch  Schmldt-Cassegraln 
receiving  telescope  end  detected  by  a  ITT  FW130  photomultiplier  operated 
ee  a  photon  counter.  The  detected  photon  bursts  are  sent  to  a  Malvern 
K7023  digital  correlator  and  a  HF85  computer  for  analysis.  The  wind 
speed  is  also  measured  by  a  cup  anemometer  with  associated  electronics. 

Tha  fen,  Indicated  in  Fig.  2,  is  used  for  the  laboratory  experiments  only. 
For  our  Initial  LDV  field  experiment  [6],  the  focussing  telescope  was 
removed  and  the  two  beams  were  made  to  cross  at  the  sensing  volume 
forming  intensity  interference  fringes.  For  field  experiments  to  verify 
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the  scaling  law,  the  beam  splitter  was  removed  and  one  laser  beam  was 
focussed  to  a  spot  In  the  sensing  volume.  The  setup  for  Indoor  laboratory 
experiments  [9,10]  is  similar  except  that  the  focussing  and  receiving 
telescopes  were  replaced  by  simple  lenses,  and  the  wind  was  created  by  a 
fan  running  at  a  constant  spaed. 

Before  conducting  serious  studies  of  laser  velocimetry,  the  above 

setup  was  used  to  clarify  the  sensitivity  of  the  photon-burst  method, 

and  determine  the  minimum  detectable  burst  counts.  For  this  purpose, 

many  LTV  and  LDV  speed  measurements  under  minimal  stray  light  conditions 

were  made  in  the  laboratory.  Six  selected  resulting  correlation  functions 

are  shown  in  Fig.  3.  Figure  3(a)  corresponds  to  a  particle  just  big 

enough  for  a  LTV  speed  measurement.  Although  the  slgnal-to-noisa  is 

not  great,  two  humps  are  clearly  discardable.  The  total  signal  count 

of  the  photon  burst  is  seen  to  be  less  than  10  counts.  Figure  3(b) 

represents  background  noise  which  has  been  somewhat  reduced  by  the 

preset  clipping  level.  Figure  3(c)  depicts  a  beautiful  LTV  correlation 

function  from  which  particle  speed  can  be  detarmined  accurately.  Even 

for  this  high  level  of  slgnal-to-noise,  only  130  signal  counts  are 

needed.  Notice  the  total  signal  plus  noise  counts  are  nearly  the  same 

for  all  three  figures  while  the  total  signal  counts  are  quite  different. 

This  suggests  chet  tho  single-particle  photon  bursts  occur  during  s 

-A 

very  short  time,  about,  A  *  10  esc.  in  tnesa  cases.  On  the  other  hand, 
the  clipped  count  gives  a  fairly  good  idea  of  the  signal  bursts. 

In  tha  absence  of  noise,  a  photon  burst  of  two  counts,  one  scat¬ 
tered  from  each  beam,  is  in  principle  enough  for  a  LTV  speed  measure¬ 
ment.  However,  even  with  clipping,  some  noise  is  present;  ten  signal 
counts  is  seen  to  be  more  than  enough  for  a  speed  measurement  in  practice. 
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The  sensitivity  of  photon  burst  correlation  is  clearly  demonstrated.  As 
•  side  Issue,  we  point  out  that  the  second  peak  in  Fig.  3(c)  Is  seen  to 
be  taller  than  half  of  the  first  peak.  This  is  due  to  the  distortion 
caused  by  the  process  of  clipping  [17].  Although  the  shape  of  the 
correlation  function  nay  be  distorted,  the  peak  separation,  and  thus  the 
measured  speed,  will  not  be  altered  by  the  process  of  single-clip 
correlation. 

Figures  3(d)  -  3(f)  are  correlation  functions  corresponding, 
respectively,  to  minimal  detectable,  noise,  and  good  LDV  speed  measure¬ 
ments.  Similar  comments  can  be  made  for  these  figures.  In  general, 
the  sensitivity  of  both  methods  are  comparable;  the  measurement  time 
used  was  the  same  for  both  experiments  although  the  laser  power  for  LDV 
measurements  was  a  factor  of  two  higher.  Comparing  Figs.  3(a),  (b),  (c) 
to  Figs.  3(d),  (e),  (f),  one  observes  without  much  surprise  that  the 
minimal  clip  count  needed  for  a  speed  measurement  is  higher  for  LDV. 
However,  once  exceeding  this  minimal  clip  count,  the  sensitivity  of 
single-particle  LDV  quickly  improves  and  catches  up  with  that  of  LTV. 

Using  the  probability  density  mode  of  the  digital  correlator,  the 
count  distribution  of  photon  bursts  can  be  measured  [9].  To  demonstrate 
the  usefulness  of  the  count  density  distribution,  we  compared  the  photon 
bursts  in  a  preset  time  interval  resulting  from  scattering  from  a 
blackened  paper  to  that  from  indoor  natural  aerosol;  the  results  are 
shown  in  Figs.  4(a)  and  4(b),  respectively.  The  distribution  of  Fig. 
4(a)  is  Poisson  indicating  the  constant  intensity  of  the  sensing  volume. 
The  gradual  tailing  off  of  the  distribution  in  Fig.  4(b)  indicates  the 
availability  of  large  single  particles,  although  with  small  probability. 


12 


IV.  The  Scaling  Law  and  Performance  Evaluation 

The  need  for  a  figure  of  merit  which  takes  Into  account  the  aerosol 

condition  to  evaluate  the  performance  of  laser  veloclmetry  has  been 

recognised,  and  the  measurement  rate  was  suggested  as  a  logical  choice 

in  our  early  work  [6],  With  the  photon-burst  correlation  technique,  we 

have  Initiated  19]  and  more  recently  succeeded  [10]  in  the  formulation 

of  a  scaling  law  for  laser  veloclmetry  in  visible  wavelengths,  relating 

tha  measurement  rate  of  atmospheric  crosswind  due  to  individual  partldaa, 

tha  signal  partial#  arrival  rata,  F,  to  laser  power  PQ  and  maaaurament 

range,  R,  under  a  glvan  aerosol  condition.  Although  two  laser  beams  are 

normally  uaad  in  laser  veloclmetry,  we  used  e  single  Gauseian-profllad 
-2 

beam  with  in  i  radius,  p,  for  tha  discussion  of  the  scaling  law.  This 
la  parmiaaibla  lines  In  a  LDV  aatup,  two  baama  are  croaied  Into  ona  spot 
containing  lntanalty  lnuarferanca  fringaa,  and  In  a  LTV  aatup,  moat 
individual  partlclas  travoraa  both  lasar  baama  because  the  beam  aspira¬ 
tion  (typically  'vcm)  is  usually  smaller  than  the  characteristic  length 
describing  the  turbulence  of  atmospheric  wind  field  [10].  As  a  result, 
a  single-beam  arrangament  should  yield  the  same  (or  nearly  tha  same) 
signal  particle  arrival  rate  (measurement  rate),  F.  Using  tha  concept 
of  minimum  detectable  photon  buret,  n^,  and  the  associated  power  depan¬ 
dent  minimum  detectable  particle  radius,  r§,  along  with  a  generalised 
Jungs  distribution  for  the  size  distribution,  dN/dr  -  cor"0,  the  scaling 
law  [10]  for  single  particle  arrival  rata  can  be  rawrittan  as 

.  ,  ^c°  t  a-B) 

(0-1)1'5  * 


F 

Avp 


(2) 
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where 


(n§/Y) 


% 


and  K0,  A,  hv,  v,  n,  i,  are,  respectively,  backucattering  coefficient 
2 

[6]  (o  »  KQr  ),  receiving  telescope  area,  photon  energy,  crosswind 
speed,  the  efficiency  of  the  detaction  system,  and  the  depth  of  view. 

If  during  a  measurement  time,  the  correlator  is  sat  to  accept  photon 
bursts  with  photon  oounts  exceeding  the  clip  level,  k,  at  least  m  times, 
the  minimum  detectable  photon  burst  [10]  is 


where  is  the  sample  time  for  the  correlation.  Typically,  m  is  set  to 
be  the  number  of  desired  observable  fringes  in  the  sensing  volume  for 
LDV,  and  mT0  approximates  the  width  of  the  photon  burst  for  LTV.  It 
has  also  bean  shown  mathematically  [10]  that  using  the  probability  mode 
of  the  digital  correlator,  the  count  distribution  of  photon  burst,  N(n), 
can  be  measured  in  a  measurement  time  Tp,  and  that  the  signal  particle 
arrival  rate  can  also  be  obtained  independently  this  way  as, 


F  -  ~  /  N(n)dn 
P  n, 

where 


(4) 


N(n) 


g  c.tvpIp  .n»-(B+l)/2 
2<S-1)0-!Y  y 
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The  correctness  of  the  scaling  law  has  been  demonstrated  in  a  laboratory 
experiment  in  which  both  signal  particle  arrival  rate  F  and  count 
distribution  N(n)  have  been  measured  independently.  The  results  are 
shown  in  Figs.  5  and  6.  Both  experiments  give  rise  to  the  same  value  of 
0  within  the  experimental  error.  The  quantity  F/(Jtvp)  in  (2)  can  be 
interpreted  as  the  number  of  particles  per  unit  volume  having  radius 
greater  than  r ,  a  quantity  familiar  to  aerosol  scientists  [19].  When 

experimental  conditions  are  such  that  a  small  range  of  r  is  encountered, 

0 

the  aerosol  parameters  cQ  and  0  are  constant  and  (2)  may  be  used  to 
determine  these  parameters  experimentally.  This  was  dona  previously  in 
the  laboratory  {10 ]  and  together  with  the  measured  photon-count  distri¬ 
bution,  our  scaling  theory  for  performance  evaluation  of  a  laser  time- 
of-f light  valocimeter  (LTV)  was  verified  in  a  self-consistent  manner. 

However,  when  a  large  range  of  rg  is  encountered,  cQ  and  0  may  not 
be  constant.  In  general,  the  measured  value  of  0  will  increase  as  the 
aerosol  radius  increases.  To  obtain  a  wide  range  of  rg  in  our  correla¬ 
tion  experiments,  the  laser  power  PQ  and  beam  radius  p  must  be  varied 
sufficiently.  In  this  manner,  the  validity  of  (2)  can  be  verified  for  a 
much  wider  range  of  rg. 

Following  the  same  procedure  outlined  in  Ref.  10  with  a  laser  power 

Po  ranging  from  0.3  raW  to  0.4.W  and  beam  radius  p  ranging  from  40  to  900 

urn  by  using  different  focussing  lenses  with  R  ■  1.7  m  and  v  ■  5.2  m/sec, 

we  have  conducted  extensive  laboratory  experiments  covering  r  values 

0 

from  0.1  to  4.0  urn.  The  experimental  results  when  plotted  as  F/(Avp) 

vs  r  is  shown  in  Fig.  7.  The  resemblance  of  Fig.  7  to  the  typical 
0 

cumulative  number  distribution  of  aerosols  [19]  is  quite  obvious.  The 
0  value  changes  from  3.2  to  6.1  as  r  increases  in  general  agreement 

9 


with  the  expected  room  aerosol  distribution.  Our  previous  result  of  ft 


4.3  occurs  gt  r(  9  1.5  im  as  a  special  case. 


Vi 
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V.  The  Effect  of  Atmospheric  Turbulence  and  Performance  Prediction 

In  order  to  use  (2)  to  predict  the  veloolmetry  performance,  a 
realistic  acaling  law  must  first  be  established  in  the  intermediate 
ranges  from  field  experlmente.  the  atmospheric  aerosol  condition  is 
different,  end  the  values  cQ  and  a  are  expected  to  be  different,  from 
ehoee  obtained  in  the  laboratory.  Atmospheric  turbulence  should  be 
monitored  and  its  affect  on  the  focuaaed  beam  radius,  p,  assessed.  With 
the  scaling  lav  tested  by  the  field  experiments,  the  performance  of  the 
laser  veioclmeter  at  long  ranges,  such  as  500  m  can  than  be  predicted. 

It  la  vail  known  that  atmospheric  turbulence,  described  by  a 
constant  denoting  random  refractive  index  fluctuations,  C*,  produces 
undesirable  effects  of  beam  wander  and  beam  breathing  [20].  Thia 
reduces  the  visibility  of  intensity  intarfaranca  fringes  in  LDV  opera¬ 
tion  [21],  In  the  summer  of  1975  during  our  LDV  field  experiments  [6], 
va  also  observed  visually  that  although  the  Interference  fringes  at  a 
range  of  1000  m,  formed  by  unfocussed  argon-ion  laser  beams,  broke  into 
patches,  clear  interference  fringas  could  be  seen  within  the  patches 
for  e  time  long  enough  for  e  photon-burst  correlation  measurement.  The 
quantitative  effect  of  turbulence  on  beam  radius  is  still  an  open 
question.  To  our  knowledge,  among  the  vast  amount  of  publications  on 
atmoapharlc  turbulence,  only  a  few  articles  have  addressed  thia  queation 
of  beam  radius  directly.  Fried  [22]  introduced  the  idea  of  optical 
raaolution  for  imaging  of  a  uniform  circular  beam  through  the  atmosphere. 
His  short  exposure  phase  structure  functions  have  been  modified  by 
Lutomirskl  et  al  [23],  and  Richardson  [24]  has  used  Fried's  optical 
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resolution  function  to  calculate  beam  radius.  We  have  recently  modified 
Fried/Richardson' h  formula  for  a  Gaussian  beam  profile  and  mudc  focussed 

beam  radius  measurements  up  to  a  range  of  100  m  using  a  photo-diode 
array.  The  experimental  radii  hava  been  compared  against  three  theo¬ 
retical  models  of  turbulence  Induced  beam  radius t  Fried/Richardson’ s 
uniform  circular  beam  [22,24]  and  the  Gaussian  beam  models  of  Lutomlrskl 
at  al  [23]  and  of  Kallay  at  al  [25].  The  detailed  analysis  and  comparisons 
will  be  published  elsewhere  [25],  but  the  main  result  is  shown  in  Fig.  8. 

It  Is  seen  that  the  model  of  Kelley  et  al  agrees  with  the  experimental 

results  better.  The  calculated  beam  radius  at  500  m  range  for  C2  -  io"^ 

n 

-2/3 

m  ,  are  -  12,340  pm,  rTC  -  12,860  pm,  rTA  »  3,572  pm,  for  the 
models  of  Fried/Richardson,  Lutomlrskl  et  al,  and  Kelley  et  al,  respec¬ 
tively.  Again,  the  model  of  Kelley  et  al  seems  to  tgree  with  the 
observed  focussed  beam  radius  at  500  m,  which  when  observable  is  about 
0.5  cm,  under  the  experimental  conditions. 

To  establish  the  scaling  law  for  field  experiments,  a  eat  of  single- 
beam  exparlmants  up  to  a  range  of  150  m  have  been  made  in  the  nights  of 
summer  1982  and  the  signal  particle  arrival  rates  monitored  by  a  computer. 
The  atmospheric  turbulence  was  simultaneously  monitored  by  a  N0AA  C* 

meter  [26].  The  beam  radii  for  the  experimental  conditions  including 
2 

measured  CR  ware  calculated  with  the  three  models  described  above  along 

with  a  Gaussian  diffraction  model  ignoring  turbulence.  The  signal 

particle  arrival  rate  F  and  minimum  detectabla  particle  radius,  r  ,  were 

8 

determined  from  the  experiment  [10,25]  and  the  data  using  the  three  sets 
of  calculated  p  are  fitted  to  (2).  This  procedure  yields  the  best  fitted 
aerosol  parameters  cQ  and  6,  shown  in  Table  I.  Although  the  different 
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turbulence  models  predict  quite  different  beam  radii  at  high  turbulence 
(see  Fig.  8),  their  effects  on  the  slope  of  the  scaling  law  determined 
by  8  Is  seen  to  be  minimal.  This  is  due  partly  to  the  self-compensating 
nature  of  the  turbulence  effect  on  the  signal  strength  (the  total  photon 
burst  counts).  The  turbulence  increases  the  beam  radius  at  the  sensing 
volume  which  decreases  the  laser  Intensity,  but  Increases  the  transit 
time  of  the  particle  across  the  laser  beam, 

Over  a  period  of  two  years,  many  sets  of  intormedlata  range  field 

experiments  Were  made  [27,28]  for  the  investigation  of  the  scaling  lav, 

but  only  the  experimental  results  of  the  summer  1982  were  used  for  the 

determination  of  cQ  and  0,  presented  in  Table  I,  because  the  values  of 
2 

Cn  were  simultaneously  recorded  for  these  data.  In  addition,  some 
single-beam  measurements  were  made  in  the  daytime  up  to  a  range  of  100 
m.  For  the  daytime  measurements ,  a  narrow-band  color  filter  (pass  band 
'vlO  A)  was  used  In  the  recaiver  to  reduce  background  light,  Due  to  the 
background,  the  signal  particle  arrival  rates  are  lees  then  the  rates  at 
night  for  comparable  ranges.  At  23  m  the  arrival  rate  was  found  to  be 
five- to- ten  times  lower. 

In  the  summer  and  fall  of  1982,  night  experiments  were  conducted  at 

Fort  Collins  to  measure  the  signal  particle  arrival  rate  at  a  range  of 
2 

500  m  with  Cn  simultaneously  monitored.  Both  single  beam  and  dual  beam 
experiments  ware  conducted  with  about  0.25  W  of  laser  power  in  the 
sensing  volume.  When  two  laser  beams  are  us? 1  the  croeswind  speed  le 
determined  by  the  ratio  of  beam  r-  deration  to  ttu  measured  time-of- 
fllght.  The  effects  of  turbulence  on  the  beam  separation  were  investi¬ 
gated  first  by  taking  a  home  movie  of  two  focussed  laser  spots  at  a 
range  of  500  m.  The  fluctuations  in  beam  separation  were  on  the  order 
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of  10Z.  Later,  the  beam  separations  were  measured  in  the  daytime  for 
ranges  of  100  m,  200  m,  and  500  m  in  the  atmosphere  using  a  linear 
diode  array.  The  fluctuation  in  beam  separation  was  typically  5Z  for 
100  m  and  200  m  range.  At  the  500  m  range,  the  measured  fluctuations 
of  beam  separation  were  from  2%  to  15%,  depending  on  the  strength  of 
the  atmospheric  turbulence. 

The  measured  signal  particle  arrival  rates  F  at  500  m  are  shown  in 

m 

Table  II  along  with  the  predicted  signal  particle  arrival  rates  and 
calculated  minimum  detectable  particle  radius  based  on  the  various 
turbulence  models. 

The  usefulness  of  including  turbulence  in  the  calculation  of  signal 

particle  arrival  rate  F  and  p  is  made  deer  by  considering  the  first  two 

experiments  (7-20-1  and  7-20-2)  of  Table  II.  These  ere  single  beam 

experiments  conducted  a  faw  days  after  the  scaling  law  experlmants  in 

the  summer  of  1982.  Predicted  F  for  the  Gaussian  diffraction  modal  is 

mors  than  an  order  of  magnitude  below  the  measured  rate  F  ,  while 

m 

predicted  F  for  the  model  of  Kelley  at.  el.  la  comparable  to  F  for 

m 

these  experiments.  The  predicted  F  for  the  other  two  turbulence  models 

is  comparable  to  F^,  but  three  to  five  times  the  predicted  F  for  the 

model  of  Kelley  et.  al.  The  remaining  two  experiments  conducted  In  the 

fall  of  1982  uaed  two  beams.  Predicted  F,  based  on  3  and  c  from  the 

0 

summer  1982  seeling  law  experiments,  is  comparable  to  F  for  all  the 

m 

models  considered.  What  is  perhaps  more  important  is  that  the  model  of 
Kelley  et.  el.  gives  the  most  reasonable  minimum  detectable  particle 
radius  among  all  models  considered.  Considering  the  Infrequency  of  the 
signal  events  (100-600  sec.  per  acceptance)  at  the  range  of  500  m,  the 
agreement  between  theory  and  experiment  is  quite  good. 
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VI.  The  Prospect  of  Routine  l.ong-Range  Operation 

As  mentioned  In  the  Introduction,  it  is  the  general  conclusion  of 
researchers  [2,5,8,16]  that  long-range  (500-1000  w)  wind  speed  measure¬ 
ments  using  cw  visible  lasers  should  be  possible  on  a  routine  basis. 
Nonetheless,  routine  operations  at  these  ranges  have  not  yet  been 
realised.  Although  the  measurements  at  500  m  reported  here  have  clearly 
demonstrated  the  possibility  and  feasibility  of  wind  speed  measurements 
in  the  visible  at  long  ranges  with  cw  lasers,  the  signal  particle 
arrival  (measurement)  rate  is  still  too  low  to  be  of  practical  and 
routine  use. 

To  further  Improve  the  performance  and  Increase  the  measurement 
rate,  higher  laser  power  should  be  used.  In  addition,  different  methods 
for  signal  analysis  which  can  handle  photon  bursts  resulting  from 
scattering  of  multiple  particles  should  be  investigated  as  well. 

Although  the  theoretical  and  experimental  measurement  rates  at  the  range 
of  500  m  were  comparable,  the  calculated  minimum  detectable  signal 
particle  radius  r  (v  200  pm) ,  is  too  large  to  be  abundant  enough  in 
normal  clear  air .  The  accepted  measurements  at  500  m  may  be  partially 
due  to  several  particles  moving  across  the  laser  beam  together.  A 
typical  correlation  function  probably  resulting  from  more  than  one 
particle  at  500  m  under  clear  air  conditions  is  shown  in  Fig.  9(a). 

This  is  compared  to  a  correlation  function  resulting  presumably  from  a 
single  particle  crossing  the  laser  beamB  at  the  same  range  in  a  drlzsle 
with  higher  wind  as  shown  in  Fig.  9(b).  The  presence  of  multiple- 
particle  signals  at  long  ranges  is  unavoidable  and  the  measurement  rate 
can  be  greatly  Increased,  if  these  signals  can  be  accepted  in.  a.  recog¬ 
nizable  form.  In  this  respect,  the  LDV  measurement  using  a  fast  Fourier 
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transform  method  [16]  which  analyzes  the  frequency  instead  of  the 
Intensity  of  a  photon  burst,  may  be  more  suitable.  In  fact,  Durst  and 
Rlchtar  have  recently  been  able  [29]  to  conduct  routine  wind-speed 
measurements  at  a  range  of  350  m.  Proper  processing  of  LTV  correlation 
functions  as  a  part  of  the  acceptance  criteria  using  some  curve-fitting 
procedure  to  enhance  the  double  humps  should  also  facilitate  LTV  meas¬ 
urements  at  long  ranges. 

In  summary,  this  report  reviews  the  photon-burst  correlation  tech¬ 
niques  of  LDV  and  LTV  for  atmospheric  crosswind  measurements  using  cw 
lasers  at  visible  wavelengths.  Measurements  based  on  photon  bursts 
resulting  from  individual  signal  particles  can  be  readily  made  up  to  a 
range  of  a  few  hundred  meters;  their  measurement  rates  can  be  understood 
by  a  theory,  including  the  effect  of  turbulence,  presented  herein. 
Experimental  evidences  of  successful  crosswind  measurements  up  to  500  m 
suggest  that  signal  resulting  from  multiple  particles  at  long  ranges 
can  no  longer  be  Ignored.  In  addition  to  the  ubb  of  higher  power  lasers 
[16,30],  proper  use  of  this  type  of  signal  and  Its  inclusion  in  the  data 
analysis  will  Increase  the  measurement  rate  at  long  ranges  (up  to  1000 
m)  and  should  make  laser  velocimetry  good  enough  for  routine  atmospheric 
applications. 
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Table  1.  Aerosol  distribution  constants  c0  and  6  as  dstsrmined  from  the 
scaling  law. 


Model  c  B 

O  !  H 


Fried/Richardson*^ 

2.83  x  10"8 

3.19  ±  0.17 

Lutomirski  et  alb^ 

1.60  x  10"8 

3.25  ±  0.18 

Kelley  et  alc) 

5.04  x  10‘10 

3.47  t  0.20 

Gaussian  Diffraction8^ 

1.48  x  IQ"10 

3.55  ±  0.21 

From  Ref.  25,  *^The  uniform  circular  model,  rT]),  b^The  Gaussian  model 
of  Lutomirski  et  al,  rTC,  c^The  modified  Gaussian  model,  rTA,  and 
d^The  diffraction  model,  r„. 


Table  II.  Experimental  results  at  a  500  ■  range. 


Fig.  1.  The  intercity  n( c )  resulting  from  the  scattering 

of  (a)  many  small  particles,  (b)  one  single  particle, 
(c)  five  signal  particles  simultaneously  traversing 
the  sensing  volume,  and  (d)  20  randomly  distributed 
signal  particles. 


Fig..  2.  Expermental  arrangement  for  LTV  measurements.  1.  Laser, 

2.  beam  splitter,  3.  focussing  telescope,  4.  fan,  5.  cup 
anemometer,  6.  anemometer  electronics,  7.  receiving 
telescope,  8.  photomultiplier  assembly,  9.  digital  rate 
meter,  10.  digital  correlator,  and  11.  computer.  The 
fan  is  used  only  in  the  laboratory  experiments  to  generate 
the  wind.  ' 


Fig.  4.  Probability  density  of  burst  counts  from  a  uniformly  illuminated 
target  (a)  and  Indoor  natural  uurosol  (b). 
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Fin.  5.  Soiling  law  of  tha  laboratory 
corralation  axparimant,  Tha 
LTV  meaauremont  rate  F  la 
plotted  as  a  function  of 
experimental  parameters, 
such  as  lager  power  P0  and 
detector  range  R.  The 
elope  of  this  plot  yields 
the  power  of  the  aerosol- 
size  distribution,  B  ■ 

4.332  ±  0.258. 


Fig.  6.  Photon-count  distribution 
N(n)  of  signal  particles 
plotted  as  a  function  of 
photon  counts  per  eample  n 
divided  by  y,  which  contains 
experimental  parameters. 

The  slope  of  the  plot  yields 
the  power  of  the  aerosol 
size  distribution,  0  ■  4.54  ± 
0.35. 
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1(mw) 


Fig.  7.  Scaling  law  of  the  laboratory  experiments  over  a  wide  range 
of  rg.  The  LTV  measurement  rate  F  divided  by  (Avp)  Is 
plotted  as  a  function  of  minimum  detectable  radius  rB  which 
depends  on  experimental  parameters  such  as  laser  power  P0 
and  detection  efficiency  n*  With  the  wide  range  of  rs 
(0.1-4  pm)  covered,  the  elope  of  this  plot  yields  varied 
values  for  the  power  of  the  aeroeol-siase  distribution,  (t  * 
3. 2-6.1. 


Fig.  8.  Measured  c~2  beam  radius, 

rm,  vs  predicted  beam  radii, 
r>j,  for  three  turbulence 
inodela. 


Fig.  9.  Measured  photon -buret 

correlation  functions  at 
500  m  under  different  atmos¬ 
pheric  conditions:  (a)  in 
clear  air  due  possibly  to 
several  particles,  and  (b) 
in  drizzle  due  possibly  to 
a  single  light  rain  drop, 
scattering  from  the  two 
laser  beams.  The  sample  time 
T0  (per  channel)  is  1.0  msec. 
The  beam  separation  is  1.9 
cm  and  3.4  cm  for  (a)  and 
(b),  respectively. 


